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Abstract This is the third article in a series of three papers on the resonance energy levels of 
anharmonic oscillators. Whereas the first two papers mainly dealt with double-well potentials 
and modifications thereof [see J. Zinn-Justin and U. D. Jentschura, Ann. Phys. (N.Y.) 313 
(2004), pp. 197 and 269], we here focus on simple even and odd anharmonic oscillators for 
arbitrary magnitude and complex phase of the coupling parameter. A unification is achieved by 
the use of 'PT-symmetry inspired dispersion relations and generalized quantization conditions 
that include instanton configurations. Higher-order formulas are provided for the oscillators of 
degrees 3 to 8, which lead to subleading corrections to the leading factorial growth of the 
perturbative coefficients describing the resonance energies. Numerical results are provided, and 
higher-order terms are found to be numerically significant. The resonances are described by 
generalized expansions involving intertwined non-analytic exponentials, logarithmic terms and 
power series. Finally, we summarize spectral properties and dispersion relations of anharmonic 
oscillators, and their interconnections. The purpose is to look at one of the classic problems 
of quantum theory from a new perspective, through which we gain systematic access to the 
phenomenologically significant higher-order terms. 
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Chapter 1 

Introduction 



We here continue our investigations on anharmonic oscillators, shifting the focus away from double-well 
and related potentials, which had been the subject of the first two papers [1] [2] of this series. Following 
ideas outlined previously in a brief exposition [3], we here elaborate on a unified theory of even and odd 
anharmonic oscillators, where the resonance energies are determined by a harmonic oscillator Hamiltonian 
which is perturbed by an even or odd monomial of the coordinate. 

Indeed, anharmonic oscillators represent a very well developed field of study, and the scientific literature on 
the subject is extensive (for a few exemplary literature references, see [H |5] [6] [71 [8] [9]). Still, it is quite 
surprising that a few basic questions regarding the anharmonic, one-dimensional oscillators have apparently 
not yet been fully addressed to the best of our knowledge. Some of these unanswered questions concern 
the general structure of resonance energies of the anharmonic oscillators; the imaginary part of a resonance 
energy is nonvanishing and describes the tunneling of the bound particle through the potential well. We recall 
that resonance energies with a nonvanishing imaginary part are relevant, e.g., for even anharmonic oscillators 
when the coupling parameter g in the interaction gx^ becomes negative {N even), and for odd anharmonic 
oscillators when the coupling parameter g in the interaction ^Jgx^^ is positive (Af is assumed to be an odd 
integer, the motivation for the square root is explained below). Here, we attempt to unify the treatment of 
the resonance eigenvalues of anharmonic oscillators of even and odd degree, with a special emphasis on the 
dispersion relations fulfilled by the resonance energies. 

We also investigate higher-order formulas ( "hyperasymptotics" ) for the oscillators of arbitrary degree. These 
have the form of generalized, nonanalytic expansions, where a nonanalytic (in the coupling constant) prefactor 
is multiplied by logarithmic and power series (this constitutes a triple expansion). Note that generalized non- 
analytic expansions, which contain expressions like exp(— A/^*") with positive constants A and 6, logarithmic 
terms like ln((7), and powers of 5, have recently been shown to be relevant to the description of energy levels 
of a number of one-dimensional quantum systems [T]|2], including problematic cases where the ground-state 
energy vanishes to all orders of perturbation theory, but the true ground-state energy is manifestly different 
from zero and dominated by a nonanalytic factor of the form exp[— l/(3.g)] (see Refs. (TU] |Tl] [121 [H])- 
Anharmonic oscillators, as shown here, also imply triple expansions for a complete description of the reso- 
nance energies. Because of their phenomenological importance, the anharmonic oscillators occupy a special 
position within quantum theory. 

Let us therefore summarize some of the most prominent questions investigated in the current report: How 
can we write down a satisfactory expansion, in powers of the coupling, describing the unstable cases where the 
coupling parameter multiplying the anharmonic perturbation attains a value that makes the levels unstable? 
Do resonance eigenvalues persist when we increase the modulus of the coupling parameter to infinity? What 
is the general behaviour of large-order perturbation theory, for arbitrary resonance energies of odd anharmonic 
oscillators of arbitrary degree? 
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Although the above questions have not yet been addressed in the literature to the best of our knowledge, 
many physical aspects of anharmonic oscillators are indeed well known. E.g., the spectrum of an even 
anharmonic oscillator (with a perturbation of the form gx^ with g > and even integer A^) is bounded from 
below and consists of discrete energy eigenvalues, when the Hamiltonian is endowed with L^(R) boundary 
conditions. The divergent perturbation series describing the energy eigenvalues of even oscillators are known 
to be Borel summable [4J. In addition, one may derive effective semi-analytic approximations to these energies 
such as renormalized strong-coupling expansions [8_, 9J in the "stable" sector of positive coupling parameter 
5 > 0, which are applicable to both weak and strong coupling. However, when g is continued into the 
complex plane and the boundary conditions are smoothly deformed in the sense of complex scaling [5j, 
then, as shown in Ref. ^6j, the resonance energy eigenvalues become complex. Indeed, the smooth functions 
En{N,g) describing the nth resonance energy eigenvalue of an even oscillator of degree N are analytic 
functions of the (complex) coupling g for — tt < argg < vr. There are branch cuts at negative g, and the 
resonance energies are of the form ReEn{N,g) ± ilmEn{N, g), where the two different signs are attained 
on different sides of the branch cut along the negative axis. 

For odd oscillators, by contrast, the situation is not as straightforward and has been discussed in the literature 
quite recently [14]. The most straightforward choice for the interaction term, namely gx^^ , does not lead 
to a satisfactory description because the transformation g ^ —g \n this case is equivalent to a parity 
transformation which changes the Hamiltonian but leaves its spectrum invariant. It may seem counterintuitive 
that the most natural definition of the perturbative term in this case involves the structure ^x^^ . In this 
case, employing the concept of PT-symmetry [TH] [161 EZl [HI [IS], one can show that the spectrum for 
£^(]R)-boundary conditions is real, discrete and bounded from below for negative coupling g < 0. Note that 
the case of negative coupling corresponds to a PT-symmetric interaction term ±i \f\g\x^'^ , where again the 
± sign of the perturbative term does not matter as far as the spectral properties are concerned. Starting 
from the case of negative g, one can then smoothly deform the boundary conditions when varying g in 
the complex plane, so that a branch cut for the energy emerges for positive coupling. Resonances (with 
a negative imaginary part of the energy) or antiresonances (with a positive imaginary part of the energy) 
appear infinitesimally displaced above and below the positive real g-axis, respectively. 

It has been known for a long time that the perturbation series describing the energy levels of anharmonic 
oscillators constitute asymptotic, divergent series. In 1971 (see Ref. [20J), the general leading asymptotics 
for the perturbative coefficients describing the energy of an arbitrary excited level of an even oscillator of 
arbitrary degree were derived. Here, we generalize this result to the case of odd oscillators of arbitrary 
degree, and for an arbitrary excited level. This generalization is made possible by a combination of (i) VT- 
symmetry-based dispersion relations derived in Ref. [14], and (ii) by generalized quantization condition that 
are conjectured in the current work. 

While the unification of even and odd anharmonic oscillators necessarily cannot lead to completely analogous 
formulas for both cases, our conjectured unified quantization conditions cover a number of cases simultane- 
ously: (i) real and positive coupling for even and the T'T-symmetric case for odd oscillators, and (ii) negative 
real coupling for even oscillators and positive real coupling for odd anharmonic oscillators. A "master table" 
with some currently known, important properties of oscillators of both parities is also provided (see below in 
Sec. 15. ip . We also emphasize the subleading corrections to various physical quantities of interest, such as 
decay rates of unstable resonances. These are provided for the anharmonic oscillators of degree 3 (cubic), 
4 (quartic), 5 (quintic), 6 (sextic), 7 (septic) and 8 (octic oscillator). 

The T'T-symmetry of the odd oscillators, for negative coupling, is used here as a tool facilitating the 
formulation of the problem leading to the general large-order asymptotics of the perturbation theory char- 
acterizing the odd oscillators. Note that T'T-symmetric Hamiltonians, which represent a natural gener- 
alization of Hermitian Hamiltonians when one allows for a non-standard scalar product of Hilbert-space 
vectors [15] [21] [22] [23l [24] , have recently found applications in many areas of physics. Examples include 
bound-state scattering theory [25l[26], physical realizations of PT-symmetric potentials in optical waveg- 
uides [27] [28] [291 [30], and the description of squeezed states of atoms [SlJ, as well as astrophysical con- 
texts [32]. Recently, the Lee model |33] and the Pais-Uhlenbeck oscillator model [34] have been shown to 
be free from problematic ghost states with negative norm, provided one endows the underlying Hilbert space 
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with a non-standard scalar product that ensures the equivalence of the model with a Hermitian theory and 
makes the time evolution unitary. Here, T'T-symmetry is applied to one of the classic problems of quantum 
theory. 

The general structure of our investigations is as follows. A short digression of the basic formulas underlying 
our investigation, and a brief discussion of the derivation of the generalized quantization conditions is pro- 
vided in Sec. 13.11 The leading-order general formulas for odd oscillators are also discussed in Sec. 13.11 We 
then proceed to the discussion of higher-order formulas ( "hyperasymptotics" ) in Sec. The anharmonic 
oscillators of the third up to the eighth degree are treated in Sees. 14.21 — 14. 71 respectively. A "master table" 
displaying some properties of even versus odd anharmonic oscillators is provided in Sec. 15.11 Numerical cal- 
culations for weak (see Sec. 15. 2p and strong coupling (see Sec. 15. 3p verify the analytic results. Conclusions 
are reserved for Sec. [6] 
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Chapter 2 



From Instantons and Dispersion 
Relations to General Formulas 



2.1 Definition of the Hamiltonians 

Before we indulge into the heart of the investigations, let us first mention a few conventions that are used 
throughout this article. The degree of an oscillator is denoted by m for a perturbation proportional to g™. 
Unless stated otherwise, we replace m if to is even, and by to — Af if m is odd. The quantum 

number of the nth excited level of a one-dimensional oscillator is denoted by n = 0, 1, 2, We denote the 

Hamiltonian of an even oscillator by HN{g), where g is the coupling, 

HM[g) = -\^ + \q^+gq'' . (2.1a) 

The resonance energies of the even oscillator are E^\g). For the Hamiltonian huig) of an odd oscillator, 
we use the convention 

hM{9) = -\^ + \q' + V~9<l'' . (2.1b) 

with resonance energies ei^^\g). There are two reasons for the appearance of the square root of the coupling 
parameter in Eq. (j2.1b|) . The first reason is that within perturbation theory, the square root ensures that 
the perturbation series for the resonance energies of the odd contain only integer powers of the coupling, 
and that the first nonvanishing perturbation is of order g, not g"^. The second reason is connected to the 
analytic structure of the dispersion relations for the resonance energies, as discussed in Sec. 12.41 below. 
Roughly speaking, the square root is necessary in order to ensure a one-to-one mapping of the spectrum of 
resonances of ()2.1b|) with a particular value of g. For purely imaginary coupling y/g = with real /3, the 
spectrum of hm^g) is real by virtue of T'T-symmetry. Furthermore, the spectrum is the same for yfg = i /? 
and y/g ~ — as well as the corresponding eigenfunctions (up to a parity transformation). Both of the 
mentioned values of y/g are uniquely associated with the value g = — and it is this one-to-one mapping 
of the resonance eigenvalue to a value of g that we need in order to formulate a dispersion relation for the 
resonance energies, which is discussed below. 

2.2 Instanton Trajectories 

Instantons are nontrivial saddle points of the Euclidean action of a quantum theory, as explained in the first 
two papers of the current series [HE]- The Euclidean action is the action obtained by varying the total energy 
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Figure 2.1: Trajectories (world lines, the time axis is the ordinate) of several instantons. The 
quartic {N ~ 4) instanton configuration Xd(^) given in Eq. ()2.11b|) and shown in Fig. (a). 
The cubic {M = 3) instanton Xci{t) is given in Eq. (|2.15bp and shown in Fig. (b). The cubic 
instanton 7;ci(^) = ^Xci(0 is found for the reflected cubic potential, where the sign of the rj'^ term 
is inverted [see Eq. (j2.20p ]. Figure (d) shows the p^^ instanton which exists for the double-well 
potential, according to Eq. (|2.23p . 



(kinetic plus potential energy) of a particle along a classical path. The partition function of the anharmonic 
oscillator can be written as 

Z{f3) - Tr exp{-/3H) , (2.2) 

where H is the Hamiltonian of the system [H — Hiq{g) for an even oscillator, and H = hM{g) for an odd 
oscillator; in our notation, we suppress the dependence of the right-hand side on the coupling parameter g]. 
A path integral representation of Z{13) involves periodic classical orbits. 



Tr exp(-/?i7) = 




(2.3) 



where the integration over the end point go of the classical trajectory (which is equivalent to its starting 
point) is explicitly given. The Euclidean action S[— j/?, ^/?, q] is the integral over the kinetic and potential 
energy of the system as it evolves from Euclidean time — 1/3 to i/3. For an even oscillator of degree N, we 
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have 



^/2 



Sw[-i/3,i/3,g] = J dt 

-/3/2 

whereas for an odd oscillator of degree A/, 

0/2 

SM[-^P,i0,q]^ J dt 

-mi 



N 



(2.4) 



M 



(2.5) 



In order to calculate the perturbative expansion of, say, the ground-state energy of an anharmonic oscillator, 
one expands about the trivial saddle point of the Euclidean action, which is given by the classical trajectory 
q{i) = 0. One writes 



P/2 li/2 



/3/2 



dt / dt'q{t)M{t,t')q{t')+g / dtg(t)'^, 



-fS/2 -13/2 
13/2 13/2 



-13/2 
13/2 



dt / dt'q{t)M{tX)q{t')^s/g / dig(t)^% 



-13/2 -/J/2 



-13/2 



M{t,t') = S{t-t'){-^ + l 



The inverse of M is the propagator 



dp e'P 1 



27r p2 



= _ exp(-lt-t'l). 



(2.6a) 

(2.6b) 
(2.6c) 

(2.7) 



As described in Ref. [35j, the perturbative expansion for the ground-state energy 



Eo{g) = lim 



i ln(Trexp(-/3ff)) 



(2.8) 



involves Feynman diagrams which can be expressed in terms of integrals of the progagator (|2.7p . Because 
the exponentiated [see Eq. (j2.2p ] quadratic quantum fluctuation about the saddle point given by the M 
operator in Eq. (|2.6p is equivalent to a Gaussian integration measure, one calls the saddle point given by 
the classical path q{t) = the Gaussian saddle point. 

In view of the obvious dominance of the Gaussian saddle point of the Euclidean action for the perturbative 
expansion, one might wonder if other nontrivial saddle points exist, and if yes, what role they might play. The 
Euclidean action corresponding to the path [q{t) = 0] is zero. Instantons are classical paths which lead to a 
finite Euclidean action in the limit of infinite Euclidean time (3. They do not correspond to a physical motion 
of the particle in the potential under investigation, but rather, to a classical motion in an inverted potential. 
We remember that the Lagrange function for a classical particle is T — y, whereas in the Euclidean action 
we are integrating over T + y (T is the kinetic energy, and V is the potential). 

In order to appreciate the importance of the instanton configurations (which are depicted in Fig. 12. ip . we 
must consider a slightly more complicated situation, namely, the cut of the partition function in the complex 
plane which is present for "unstable" values of the coupling parameter g, i.e., for those which lead to a 
potential with unstable, relative minima and potential wells through which the particle may tunnel [36, 37j. 
For the even oscillators as given by the Hamiltonian (|2.1ap . the unstable coupling parameters are given 
by g < 0, because the potential ^q"^ — \g\q^ approaches — oo for q ±oo. For an odd oscillator [see 



Eq. ()2.1bp ]. all positive values of g lead to an unstable potential because ^q'^ + gq^^ — oo for q —oo. 
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The ground-state energy in this case acquires an imaginary part and becomes a resonance energy. The 
imaginary part of the resonance energy is related to the imaginary part of the partition function according 

The sign of the imaginary part of Eo depends on the sign of the (infinitesimal) imaginary part of g. Discussing 
(for definiteness) the case of an even oscillator, we define Eo{g) to be the resonance energy [with Imi?o(.g) < 
0] for 5 = — |,g| +iO and to be the anti-resonance energy [with lmEo{g) > 0] for g — ~\g\—iO. The cut of the 
resonance energy then coincides with the cut of the partition function. It has been observed in Refs. [3711351 
that the contribution of the Gaussian saddle point to the cut of the partition function, i.e., to lmZ{(3), 
vanishes, and that lmZ{/3) can be approximated very well by the instanton configurations. 

In order to find the instanton configuration for the even oscillator, we scale the classical path as q{t) — 
(-5)^1/(^^2) in Eq. (jZ?]) . The Euclidean action becomes 

(-g)-w^ /dt fU^ + U'-e^) . (2.10) 



The instanton now is to be searched in the space of the paths ^{t), whose Euclidean action remains finite 
when t goes from — oo to oo. For = 4 and to = 0, the motion corresponds to a classical trajectory in the 
inverted potential ^"^ — and is depicted in Fig. I2.1f a'). It becomes immediately clear that the instanton 
configuration is twofold degenerate. Namely, starting for Euclidean time i ^ oo at a position ^{t = —oo) 
infinitesimally to the left or to the right of the central maximum of the inverted potential, the particle 
may go through either one of the two "troughs" to the right or to left, before bouncing back toward the 
central maximum at Euclidean time t ^ oo. According to Refs. [SHlES], the corresponding solutions to the 
(Euclidean) equations of motion read, for general N, 

'ZdW = (-ff)"^X?iW, (2.11a) 

Here, to is the "starting point" of the instanton, acting as a collective coordinate for the instanton config- 
uration from the point of view of the path integral formalism \35j. The instanton action is 

2x2/(7V-2) / ^ N \ / i\2/{N-2) 



A(m) = 22/("-2) ^ f 



(2.13) 



m — 2 m — 2^ 

where B{x, y) = r{x) r{y)/r{x + y) is the Euler Beta function. 

For odd anharmonic oscillators of degree M, we transform q{t) — — g^^/^^*^^'*^x(t) and obtain the scaled 
Euclidean action [cf. Eq. (|2.5p ] 

SM=5-^ j At{\x^ + \X^-X'') ■ (2.14) 

With this scaling, the inverted potential now reads x*^ — and is of the same form as for the even 
oscillator. The instanton trajectory q{t) = qc\{t) is unique and reads 

qcx{t)=9'^Xci{t), (2.15a) 

^^'^^^ " {l + cosh[(A/f-2)(f-to)]}i/(M-2) ■ (2.15b) 



10 
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For the case M = 3 and t ~ to, the trajectory is shown in Fig. I2.1f b'). Inserting qci{t) into Eq. (|2.14p . we 
obtain the instanton action 



Both instanton actions (|2.12p and (|2.16p are positive in the domain where the instantons exist (for negative 
coupling in the case of an even oscillator and for positive coupling in the case of an odd oscillator). 

The formulas (|2.12p and ()2.16p can be unified. Namely, if we parameterize the perturbation as gq™ for 
general (even or odd) m, then we can write the instanton action as (T2/.g)^/''"~^^ i?(;;f32 ' 7^1^)' where 
the plus and minus signs correspond to odd and even potentials, respectively. 

The scaled instanton action A{m) is expressible in terms of a scaled potential W{m,q) where the sign of 
the perturbative term has been inverted, 

>V(m,g) = ig2_^™ (2.17) 

As evident from Fig. 12.11 the instanton configurations correspond to the motion of a classical particle in a 
potential —W{m,q). We obtain from the paths given in (|2.11bp and (|2.15bp . 

2-l/(m-2) 

Aim) = j dq (2 v/2W(m,,)) = 2^'^^-^^ B f ^) , (2.18) 

J V / \m — 2 in — 2 J 



where we note that 2"-^/^™"^) is a zero of the potential yv(m,(7). 

The instanton path defined in Eq. (|2.15ap fulfills qc\{t) > 0. However, this instanton configuration is not 
uniquely associated with the Hamiltonian (|2.1bp . To see this, we start from an odd anharmonic oscillators 
of degree M and transform q{t) — (^^1/(2*^-4) j^^^y |p ^f^i^ manner, we obtain the scaled Euclidean action 
[cf. Eq. (ESI)] 

SM - ff-^ Jdt{^r,'+'2V' + V"") ■ (2.19) 

With this scaling, the inverted potential reads — (77*^ + ^ 77^) = (— ?;)^^ — ^ 77^. By comparison to Eq. (|2.15bp . 
the instanton is found to be 

= " {l + coshp./-2)(^-M]}^/(^-^-^) " ■ ^'-'^^ 



For M = 3 and to = 0, the corresponding instanton trajectory is shown in Fig. I2.ir cl. 

For the double-well potential, the Euclidean action reads, according to Eq. (2.10a) of Ref. [1], 

Sdw - / + ^ - P)') ■ (2.21) 

This expression contains the double-well potential 

Vd„(<7)-ip'(l-p)^ (2.22) 



One of the possible instanton configurations is 

't-to 



ptiit) = exp 



9 



(2.23a) 



where again to is a collective coordinate. The instanton configuration joins the two minima of the double-well 
potential. It does not define a periodic path whose starting point at t — —00 is the same as the end point 
at t ~ 00. In order to construct a periodic orbit, one has to 'join" two instantons according to Fig. 5 of 
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Sec. 4.4.1 of Ref. [IJ, so that the pseudoparticle returns to its initial point. This can be done by joining the 
instanton trajectory ply{t) with the trajectory p^i{t), where 

+ l) , (2.23b) 
which is shown in Fig. I2.1f d). for the case Iq — 0. The instanton action is found to be 

Sdwk'i]^^- (2-24) 

Let us recall an important difference of the odd oscillators (and of the double well, by the way) to the 
even oscillators. In the former cases, the instanton configuration is found for all positive, real values of the 
coupling parameter. In the latter case (even oscillator), one first has to scale the action to be consistent 
with the case of negative coupling [see Eq. (|2.10p ]. before solving for the instanton [see Eq. (j2.15bp and 
Fig. I2.ir a')]. For the even oscillator, the instanton only exists for specific regions in the space of complex 
coupling constants. 



fexp 


to-t 
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2.3 Quantum Fluctuations and Decay Widths 



According to Refs. [37] [35], the quantum fluctuations about the instanton configurations give rise to a 
functional determinant of a Bargmann potential, whose spectrum and determinant are analytically calculable. 
For the ground state of an even oscillator, this calculation leads to the result 

1/2 

exp 



Im4^)(.g) 



1 / 2C{N) 

l^(_g)2/(JV-2) 



A{N) 



where 



C(m) 



1 



22/(m-2) 



exp 



dq 



V2 



(_^)2/(W-2) 



_ 22/("i-2) 



(2.25) 



(2.26) 



A slight generalization of this result to the nth resonance energy E,, 
leads to 



Imi5W(5) 



1 



2C{N) 



!V2^ V (-5)2/(^-2' 



ri+l/2 



'"^\g), for potentials of even order N, 
A{N) 



exp 



(_g)2/(iV- 



(2.27) 



We take the opportunity to correct a misprint in Eq. (8) of Ref. [3], where the sign of the mantissa of the 
exponent n + 1/2 was accidentally inverted. This result is valid in leading order; corrections are of relative 
order g'^/i^^'^) . For odd potentials in the normalization given by Eq. (|2.1bp and for positive coupling, we 
have for g ^ 0, 



Im e, 



(5) 



1 



2C{M) 



1/2 



exp 



A{M) 
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l/(Af-2) 



The generalization to arbitrary n reads 



1 



2n!V27r \9 



2C(M) 



jl/(M-2) 



n+1/2 



exp 



A{M) 

' gl/{M-2) 



(2.28) 



(2.29) 



The sign of the imaginary parts of the resonance energies in Eqs. (j2.25p . ()2.27p . ()2.28p and ()2.29p is 
negative, as it should be for a resonance energy (as opposed to an anti-resonance, where the imaginary 
part of the energy is positive). As outlined in Refs. [371 [35], the sign of the imaginary part of the energy is 
somewhat arbitrary as it is determined by the square root y/eo = ±i|y^|, where cq < is the (only) negative 
eigenvalue of the matrix that describes the quantum fluctuations about the instanton configurations. As we 



see in the following, the resonance energies E^\g) and en"'{g) actually have cuts along the negative and 
the positive real g axis, respectively. We have to associate either sign of the imaginary part with a particular 
side of the branch cut. This becomes clear when we discuss the dispersion relations fulfilled by the resonance 
energies. 
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2.4 Dispersion Relations 

The dispersion relations formulate connections of the values of the resonance energies E^\g) and ei!^^\g) 
of the even and odd oscillators as a function of g in the complex g plane. Before we discuss them, let us 
first justify the use of the convention y^q^' for the perturbative term in the odd oscillator of degree M [see 
Eq. ()2.1bp ]. The dispersion relation for odd oscillators is intimately linked with the square-root convention 
for the coupling term. 

For illustrative purposes, we consider the Hamiltonian in the normalization hM{g) = —\dq + ^'f' + y/gi^^ . 
with odd M. A change in the sign of the coupling term ^ — > —^/g can be compensated by a parity 
transformation q —q. So, if V'(q) is a resonance eigenfunction of the cubic Hamiltonian, suitably continued 
into the complex plane, then i^{—q) is a resonance eigenfunction for the cubic Hamiltonian with a sign 
change of y^, but the same resonance eigenvalue. We conclude that the resonance spectrum of the cubic 
Hamiltonian is invariant under the transformation y/g ^ —\f§- 

It is thus natural to formulate the cubic Hamiltonian as a function of because once g is specified, 
the ambiguity in defining the square root is absorbed into the invariance of the resonance spectrum of the 
Hamiltonian under a sign change of the coupling term. A specification of g alone suffices to define the 
spectrum of resonances. As a function of g, the resonance eigenvalues of have branch cuts along the 

positive real axis. For g ~ Re((7) + iO, the imaginary parts are negative, whereas for g — Re(.g) — iO, the 
imaginary parts are positive. For negative coupling g — —\g\ < 0, the resonance eigenvalues of /i3((?), 

h^{~\g\) = -\dl^\q^±\\g\^l^q\ g>0, (2.30) 

have no imaginary part at all because the cubic Hamiltonian is T'T-symmetric for imaginary coupling and has 
a purely real spectrum. Therefore the only branch cut of the resonance eigenvalues of the cubic Hamiltonian 
is along the positive real axis. Note that the T'T-symmetry of the cubic Hamiltonian here serves as a property 
which facilitates our analysis of the structure of the branch cuts. 

In general, the dispersion relation [14J for the resonance energy eigenvalues em^\g) of a general odd oscillator 
therefore reads 

'[9)^n+- + - ds ^ . (2.31a) 

2 TT J s{s-g) 



The subtracted dispersion relation [39] [20l [6] for the energies En{N,g) of the even anharmonic oscillators 
of degree N is 

Er(,)-n+l-i /d., -°5";'°' . (2.31b) 

2 TT J s{s-g) 

— oo 

These are subtracted dispersion relations; the term n + ^ fixes the energy level for vanishing coupling. 



2.5 Generalized Bender-Wu Formulas 

The general expressions for the imaginary parts of the resonance energies given in Eqs. (|2.27p and (|2.29p . 
together with the dispersion relations ()2.31ap and (|2.31bp . can be used in order to derive integral repre- 
sentations for the perturbative coefficients and to derive their leading factorial growth for large orders of 
perturbation theory. Based on Eq. (|2.31bp . we can expand the integrand in powers of g, and obtain the 
coefficients of the perturbation series for the energy levels of an even oscillator. We find for the nth level of 
the oscillator of order N, 
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2.5. GENERALIZED BENDER-WU FORMULAS 



Because Im£^i^'(s + iO) < for — oo < s < and because the integration variable s is negative along the 
entire integration contour, the sign of the perturbative coefficients is determined according to the formula 
(— 1)^+-^ for K > 1. By contrast, for the nth level of an odd anharmonic oscillator of degree M, we have 



AM) 



g(M) K 



AM) _ 
^n.K>l ~ 



K=0 



- I ds 

TT 





Im£l^'^(g + iO) 

„K+1 



(2.32b) 



We can again argue that lmel^^\s + iO) < for < s < oo, i.e. along the entire integration contour. This 
implies that the perturbative coefficients for the odd anharmonic oscillators are all negative, except for the 
zeroth-order term (in g), which is simply n + ^. 



Based on Eq. ()2.27p and (|2.31bp . it is a simple matter of integration to calculate the large-order perturbative 
expansion for an arbitrary level of an even oscillator of arbitrary degree. 



E. 



(N) 



(_l)if+i(iV-2) 



r 



-K 



B 



N 



N 



N -2' N -2 



K 



(2.33) 



which is valid up to corrections of relative order K ^. We here confirm the result of Bender and Wu [20]. 
Our new result concerns the case of odd oscillators, where we integrate (j2.31bp using Eqs. (j2.29p and find 



-n,K 



M 



7r3/2 n\ 2^^+ 



— r{{M-2)K 



B 



M 



M 



Af - 2 ' M - 2 



-{M-2)K-n- 



K OO. 

(2.34) 



The formulas for the KtU order perturbative coefficients e'^^j^ and e^*^ have been compared to numerical 
calculations for a number of example cases, as discussed below in Sec. |4l 
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Chapter 3 



From WKB Expansions to Generalized 
Quantization Conditions 



3.1 Perturbative Expansion 

The approach outlined in Sec. [2] is inspired by field-theoretical considerations and leads to an intuitive un- 
derstanding of the physics involved in the derivation of the leading-order factorial growth of the perturbative 
coefficients, as given in Eqs. (|2.33p and (|2.34p . However, this formalism does not make use of the additional 
simple structure of the problem at hand, namely, the structure of a one-dimensional nonrelativistic quantum 
mechanical oscillator for which semiclassical or Wentzel-Kramers-Brioullin (WKB) and perturbative expan- 
sions can be derived without any further technical difficulties. Here, we strive to combine both methods 
in order to derive higher-order formulas which grant us access to phenomenologically relevant, numerically 
large correction terms. 

To this end, we first investigate the derivation of higher-order perturbative approximants to the wave func- 
tions and to the energy eigenvalues of the anharmonic oscillators. Suitable scaling transformations in the 
original Schrodinger equation corresponding to the Hamiltonians ()2.1ap and ()2.1bp prove useful in that 
regard. For the even oscillator given by Eq. (j2.1ap . we start from 




(3.1) 



and transform according to q 



9 



1/{N~2) 



q, which leads to 




(3.2) 



For (j2.1b|) . we start from 




(3.3) 



and transform according to q 



9 



l/(2(Af-2)) 




(3.4) 
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3.1. PERTURBATIVE EXPANSION 



We identify tfie general coupling parameter 

G^Girn,g)=g-i-\ X (m) ^ '-t^I^l^^ , (3.5) 

TO — 2 

which implies that G{M,g) = 51/(^-2) j^j G{N,g) = g2/(A'-2) g^g^ ^ (^g^g^ jx| is the 

largest integer smaller than X). With this identification, Eqs. (|3.2p and p.4p have the general structure 

V{m,q)^^q' + q"'. (3.6) 

Note that the potential V(m, g) differs from the potential yv(m, g) used for the calculation of the instanton 
configurations, in the sign of the term [for the definition of W{m, q), see Eq. (j2.17p ]. 

Using the scaled and unified form for the Schrodinger equation given by Eq. p.6p . the wave functions and 
resonance energies are amenable to a perturbative treatment. The general formalism has been outlined in 
Sec. 3 of Ref. f2] and is recalled here, with a special emphasis on the anharmonic oscillators. We first 
transform the Schrodinger equation to the Riccati equation by setting 

^^.fiql_s^ 
^{q) G ' ^{q) G2 G ' ' 

The Riccati form of the Schrodinger equation then reads, 

Gs'iq)^ s^{q)+U\q) = 0, (3.8a) 



U{q) = v/2[V(TO,g)-G£;]. (3.8b) 

The coupling G formally takes the role of ft in the scaled equation. The function S{q) is proportional to the 
logarithmic derivative of the wave function. It has an implicit dependence on the degree to of the oscillator, 
on the coupling G and on the energy E, which we suppress in our notation. We can uniquely decompose 
S{q) into even-parity and odd-parity (in G) components, 

s{q) = s+{q) +s-{q) = s+{m,q,E,G) + S-{m,q, E,G) , (3.9a) 

where 

s±{m,q,-G,-E) = ±s±{m,q,G,E) . (3.9b) 
This leads to the following system of equations, 

Gs'_{q) - s+{qf - s^{qf+U^{q) = 0, (3.10a) 

Gs\{q)-2s+{q)s^{q) = Q. (3.10b) 

From the second of these equations, we have 

s-{q)^i^ — TT, (3.11) 

2 s+{q) 



SO that we can write s as a function of s+ only. By virtue of the decomposition (|3.9p . this means that 

1 



^{q) = -jl= exp ( -G-' I dqs+{q) ] . (3.12) 



If we approximate s{q) w s+{q) w ^2V(to, q), then this is just the textbook version of the Wentzel- 
Kramers-Brioullin (WKB) approximation [40J. The advantage of the above formalism is that it allows for 
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3.1. PERTURBATIVE EXPANSION 



systematic expansions both for fixed E in powers of G, which leads to the perturbative expansion of the 
energy levels, but also in powers of G at G E fixed, which is an expansion relevant a priori for excited states, 
because G is assumed to be small and E must be large for G E to attain values of order unity. The latter 
approach leads to the WKB expansion. 

We start, though, with the perturbative expansion. Because the wave function for the nth excited state has 
n zeros on the real line, we can formulate the following quantization condition, 

1 I dz^ = --^ i dzs{z)^n, (3.13) 
27ri Jc ^[z) 27nG Jc 

where C is a contour that encloses all the zeros in the counterclockwise (mathematically positive) direction 
and the definition in Eq. (|3.7p for the logarithmic derivative of the wave function has been used. Equation 
(|3.13p is the Bohr-Sommerfeld quantization condition. 

We now seek a solution of Eq. ()3.10p by expansion in G at _B fixed, taking advantage of the elimination 
of S-{q) via ()3.10bp . It is easy to see that the solutions for s{q) and s+(g) can be written in terms of the 
expression 

u{q) = ^2V{m,q) (3.14) 

alone. Furthermore, the zeros of the denominators of S{q) and S+{q) are given exclusively by the zero(s) of 
u{q). Because of the simple structure of the Riccati equation (|3.8p . one can write a recursive scheme, within 
the perturbative expansion in powers of g at fixed E, that leads to higher-order perturbative approximants 
which we would like to label by a lowercase s, 

s{q) = V SK{q) , So{q) = u{q) , s,{q) - ""'^^^ " 



K=0 



2u{q) 



SK[q) 



K-l \ 
^'k-M - E SK-M)si{q) , K>2. (3.15) 
1=1 / 



2u{q) 

To order G^, we find for the symmetric part sj^{q) which is invariant under the exchange g — > —g, E -> —E, 

E 2 (}_ ""(g) _ 3 u'jqf _ 1 
u{q) ^ V 4 8 u^{q) 2 v?{q)^ 

This structure implies that the only location relevant for the singularities in the contour integral p.l3p is the 
origin in the complex z-plane. For even anharmonic oscillators, this is immediately obvious because the only 
zero of the function u{q) = q + 2q^ then lies on the real axis {N even). For odd anharmonic oscillators, 
the function u{q) = q\Jl + 2q^'' has an additional zero on the real axis at xm — —2^/*^. However, the 
residue at xm vanishes. 

Because s^{q) = ^ G s'_^{q)/s+{q) and s{q) = s+iq) + s-{q), we can finally write the quantization condition 
([TT31) as 

-2^//^^-(^)^"+i ^'-''^ 

Here, G is a contour that encloses only the origin. The advantage of this formalism is that the function 
B{m,E,G). defined by 

B{m,E,G) = -—^(fdzs+{z), (3.18) 
27riG Jc 

provides for a universal means of determining the perturbative expansion for an arbitrary excited level, by 
means of the perturbative quantization condition 

B{m,E,G) =n+^. (3.19) 

The rational is to construct higher-order approximants to the wave function according to Eq. p.isp . then 
to find the residue of the higher-order approximants to the wave function at the origin. The sum of the 
residues of the higher-order (in G) terms then defines the function B{m, E, G). 
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3.2. EXAMPLE: PERTURBATIVE EXPANSION OF THE SEXTIC POTENTIAL 



3.2 Example: Perturbative Expansion of the Sextic Potential 

The generalized coupling parameter for the sextic potential with iV = 6 is given by Eq. (|3.5p as G = g^^^. 
The zeroth-order approximant is 

So{q)^u{q) = qVl + 2q^. (3.20) 
Entering into the recursive system ()3.15p and calculating the residue at g = 0, we obtain 

Bi6,E,G)=E-G' (jE+Ie' 

V 256 32 16 I ^ ^ J 

Note, in particular, that the term -E/u{q) in Eq. ([STTe]) leads to the leading term E + 0{G^) in Eq. ((T2T|) 
in view of u{q) = q + 0{q'^). Note also that the term of order G in Eq. p.l6p does not produce any 
contributions to the residue at the origin. By explicitly solving the equation i?(6, Eff\G) = ^ for the ground 
state resonance energy i?Q^^ order by order in = g, one obtains the perturbative expansion (as a function 
ofg). 

+ f G^-^G^ + 0(G^) 

1 15 3495 o 
= -^+^9~^9' + Oig^), (3.22) 

for the ground state of the sextic potential. 



3.3 Semiclassical or WKB Expansion 

We now turn our attention to the WKB expansion, which is an expansion in powers of G for fixed GE. 
A priori, a natural starting point for the WKB expansion would be Eq. (j3.8bp . with S{q) « i^{q) = 
y^2V{m,q) — 2GE. However, a second purpose of our construction of the WKB expansion is the cal- 
culation of the instanton action, and of perturbations thereof. We have already seen in Sec. [2] that the 
instanton configuration does not exist for arbitrary values of the coupling parameter. In the case of an 
even oscillator, it exists only for negative coupling g < [see Fig. I2.ir a1]. In the case of an odd oscilla- 
tor, the instanton configuration involves positive values of the coordinates only if we assume that y/g < 
[see Fig. I2.1f bl]. For y/g > 0, the instanton configuration entails negative values of the coordinate [see 
Fig. I2.1f cl]. We therefore have to scale the coupling constants differently as compared to Sec. 13.11 if we 
wish to grasp the instanton configuration within our semiclassical expansion. 

Let us start with the case of an even oscillator, with g — ~\g\ < 0. Starting from [see Eq. (|2.1ap ]. we obtain 

+ ^-Ecb (3.23) 



and transform according to q ^ \g\ ^/^^ q, which leads to 

■q' ~qA ^^\g\^2 E<l,. (3.24) 



\9\~' d I, 



2 c)g2 ■ 2 
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3.3. SEMICLASSICAL OR WKB EXPANSION 



The sign of the term is inverted as compared to Eq. (|3.2p . For an odd oscillator, we start from 



and transform according to g — > — g 2)] which gives 

We define the general coupling parameter 

G^G(m,g) = |5|^('"), X{m) = 



gM-2 Ec 



2[[m/2] -m + 2 
m - 2 



(3.25) 



(3.26) 



(3.27) 



where X{m) is given in Eq. (|3.5p . With this identification, Eqs. (|3.24p and (|3.26p have the general structure 



1 



d 



(3.28) 



Here, we reencounter the potential W{m,q) where the sign of the term is inverted as compared to 
V{m,q). We set 

£^_%) 
G ' 

where the symbol S (as compared to S) denotes the logarithmic derivative of the scaled wave function for 
the potential W{m,q). The Riccati form of the Schrodinger equation then reads, in terms of S, 



(3.29) 



GS'{q)-S'{q)+r\q)^0, 



T{q) = y2 \W{m,q)-GE 
The recursive scheme for calculating higher-order WKB approximants then reads 



(3.30a) 



(3.30b) 



S(g) ^ G^SK(g), So(g)=r(<7), Si(g) 



K=0 



2So(g) ' 



K-l 



(3.31) 



1=1 



Unlike the perturbative approximants, which have an isolated pole at g = 0, the WKB approximants have 
a branch cut, which coincides with the branch cut of the square root function of the zeroth-order WKB 
expansion. 



So (9) 



W{m,q)-GE 



(3.32) 



In practical calculations, it turns out to be convenient to define the square root function so that its branch 
cut is along the positive real axis (this aspect will be discussed below). 

Let us here recall the perturbative quantization condition derived for the perturbative approximant s+(z) as 
given above in Eqs. (|3.17p and (|3.18p . 



* dzs+{z) ^ ~2Tri B{m,E,G) = -27ri ( n + - 
G Jc V 2 



(3.33) 
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3.3. SEMICLASSICAL OR WKB EXPANSION 



Here, C was defined as a contour that encircles all poles of s+{z) on the real axis in the mathematically 
positive sense. 

These correspond to the perturbative functions at a particular value of G. In order to make the connection 
to the formalism given above, we must therefore identify the particular value of G that corresponds to the 
same value of the original coupling parameter g as the value of G we are using for our calculations. 

For even potentials, this identification proceeds as follows. We have 

G = 52/(JV-2) (3 34) 

and G, for negative g, is given by 

G=|<?p/(^-^) = (-.g)^/(^-^), 5<0, (3.35) 

so that 

G = (-G(--)/f '^"-^^ = cxp (^) G ^ i-g f^--^^ , (3.36) 
on the cut, i.e. for g < 0. For odd M, we have 

G = .gi/(^-2) (3.37) 
G=|g|i/(*^-2) =.91/(^^-2), g>0, (3.38) 

so that 

G = G = 5i/(M-2) (3 39) 

We thus have 

m even 



G = cxp[i$(m)] G, $(m) = <( m-2 • (3.40) 

1 m odd 

Consequently, it appears useful to consider the contour integral of the WKB expansion around its cut, 

i <fdzS+{z), (3.41) 
G Jc 

where C is a contour that encircles the cut of the WKB expansion in the mathematically negative, clockwise 
sense, and S+{z) is the part of the WKB approximation S{z) that is symmetric under the transformation 
g —g, E —E. Namely, in analogy to (|3.9p . we define the decomposition 

S(g) = S+{q) + S_ {q) = S+(m, q, E, G) + S_ (m, E, G) , (3.42a) 

where 

S± (m, q, -E, -G) = ± S± (m, q, E, G) . (3.42b) 

The contour integral about the cut of the function S+ leads to a more complex structure as compared to 
the right-hand side of Eq. ()3.17p . 



i j>dz?,+ {q) = A{m,E,G) + \ \n{2T:) ^\n[r {\ ^ B{m,E,G))] +B{m,E,G) In f - 



G 



2C{m) 



(3.43) 



The left-hand side of the above equation is calculated as a function of G, a priori [including the WKB 
approximation to the wave function, S(g)]. The right-hand side defines, implicitly, functions of G. The 
procedure is to calculate the integral of S+{q) about the cut, first of all, as a function of G, 

f{G) = i j^dzS+iq), (3.44) 
c 
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3.3. SEMICLASSICAL OR WKB EXPANSION 



and then to identify 

/(G) = / (a>(m) G) = A{m, E,G) + \ ln(2^) 

- In [r (i - B{m, E, G))] + B{m, E, G) In (^-2^pW^^ . (3.45) 

Let us analyze the right-hand side of Eq. (|3.45p . To this end, we expand the expression 
ln[_r(i — B{'m,E,G))\ for large principal quantum numbers. Because 

lnr(i + z) wzln(z)-z, (3.46) 

and in view of 

B{m,E,G)^E + 0{G), (3.47) 

we can approximate the logarithmic terms on the right-hand side of p.45p as follows. Taking only the first 
term from Eq. (|3.46p . we obtain 

^\nr{\-B{m,E,G))+B{m,E,G)\n' ^ ' 



2C(m) / 

« B{m, E, G) H-E) + B{m, E, G) In (yl^^^ « B{m, E, G) In [j^^ ■ (3-48) 

This expression is real if G is positive (if E is real). However, if the arguments of both logarithms in the 
second line of the above equation acquire an infinitesimal negative imaginary part, then both logarithms 
acquire an imaginary part of ±7r«, and the term 2niB{m, E,G), previously encountered the right-hand side 
of Eq. p.33p . is recovered, thus establishing a connection of the WKB and the perturbative contour integrals. 

The strategy now is to calculate the contour integral in Eq. (|3.44p and then to use the already determined 
result for the function B{m, E, G) (see Sec. l3.1l of the current work) to subtract perturbative terms from the 
right-hand side of Eq. (|3.45p . in order to obtain subsequently higher-order terms of the function A{m, E, G) 
in the WKB limit G and G E fixed. A well-defined analytic procedure for the evaluation of the contour 
integral in Eq. (|3.44p . based on the method of asymptotic matching (see Sec. 7.4 of Ref. [41]), is presented 
below in Sec. 13.41 In leading order, one can show that 



A{m) 



Ai^,E.G).^4 '^o'^Hr' . (3.49) 



m even 



G 1 A{m) 



odd 



where A{m) is the instanton action defined in Eq. (|2.18p . We will therefore call A{m,E,G) the instanton 
function in the current article. 

At the end of the calculation, it is thus natural to reexpress the "perturbative B function" and the "instanton 
A function" in terms of the original coupling g given in Eqs. (|2.1ap and (j2.1bp . We thus define 

2([m/2]]-m + 2 

A^{E,g) = Aim,E,g — ), (3.50a) 

2[[m/21]-77i + 2 

B„,{E,g)^B{m,E,g — ). (3.50b) 

In this convention, the Bm functions have an argument of integer power (just g) and contain also only 
integer powers of the coupling in their expansion in g. The Am functions, though, contain fractional powers 
except for the oscillators of the third, the fourth, and the sixth degree (the latter one is a special case). 
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3.4. EXAMPLE: CONTOUR INTEGRALS FOR THE QUARTIC POTENTIAL 



3.4 Example: Contour Integrals for the Quartic Potential 

We remember that the contour integral around the cut of the WKB expansion is to be evaluated for the case 
where instantons persist, i.e., for g <0. We have iV = 4 and therefore, in view of Eqs. (|3.5p and p.27p . 

G=-G= -g {N^4). (3.51) 

The task is to evaluate the contour integral of the symmetric part S+{q) of the WKB expansion of the 
logarithmic derivative of the wave function, as defined in Eq. ()3.4ip . We approximate 

T = i (j^ dqS+{q) ^ -- jdq [So(<?) + 5' 83(9) + / 84(9)] = Tq + + T4 . (3.52) 
c c 

Here, the WKB approximants So(g), 82(9) and 84(9) are defined in Eq. (|3.3ip . and the integrals To, T2 and 
T4 are defined in the obvious manner. One can easily show that the contributions from the terms 81(17) and 
83(9) vanish, and we recall that C is a contour that encircles the cut of the WKB expansion in the clockwise 
direction. The zeroth-order term is 

ro = --^dg8o(<?), (3.53) 
c 

with the zeroth-order WKB approximant 



80(g) = V q)~GE= yjq^ - 2q^ + gE . (3.54) 

We define the square root function to have its branch cut along the positive real axis. Directly above the 
cut, the value of the square root is positive, while directly below, it is negative. As the contour is clockwise, 
we thus have an integration interval from zero to 1/V2, while below the real axis, we go from l/\/2 to zero. 
With this definition, the contour integral is 

1 

To = <j) dq^/f~2^fTgE = --Re J dq{q^ - 2q^ + gE)^/^ . (3.55) 
c 

The specification of the real part is necessary because we evaluate the contour integral about the cut of 
of the WKB expansion. For the evaluation, we use the method of asymptotic matching (see Ch. 7.4 of 
Ref. [41]) which involves an overlapping parameter e, that is also used in Lamb shift calculations [42l[43ll44] 
in order to separate the low-energy from the high-energy contribution to the bound-electron self-energy. The 
overlapping parameter fulfills 

0<.g£;<e, (3.56) 

and we divide the integration interval (0, l/\/2) into three parts, which are the intervals (0, e), (e, ~ e). 

and (-75 — e, 1/V2). Each result is first expanded in gE, then in e, and the divergent terms in e cancel at the 
end of the calculation. We remember that the integral is evaluated for g < 0. The integral To is recovered 
as 

ro = -- (/1+/2+/3) , (3.57) 
9 
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3.4. EXAMPLE: CONTOUR INTEGRALS FOR THE QUARTIC POTENTIAL 



where the Ij comprise the above mentioned three intervals, respectively. The first contribution is, up to 
order (gE)^, 



Ii=RcJ dq {q^ - 2q'^ + gE f/'^ 



+ {gEf 
+ igEf 



2 \gE 



{gEf 



319 

"Is" ' 

23189 



1 

5 



5 

4^2 



35 



315 



— T In 

4^2 4 



2e2 

gE 

1155 



2e2 

gE 



384 48e6 326* 326^ 



32 



In 



2f 

gE 



(3.58) 



Here, we have taken into account the fact that only the real part of the result matters; hence, the argument 
of the logarithms reads — 2e^/(gi?) > 0, because 5 < 0. The second interval leads to 



1 

V2" 



Re 
1 

6^ 



dq{q'~2q^ 



^{gE) 
{gEf 
[gEf 



.ln(- 



1 
2 
319 

"Is" 

13327 
192 
1 



gE)"' 
+ {gE? 



1 

8? 
5 



5 1 

4 " 4^ 

5 1 

4? + 6 2i/4e3/2 
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The third interval gives rise to the following results, 
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/3 = Re j Aq{q' -2q^+gEY'' ^{gEf 
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(3.59) 
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8 23/4g5/2 96 2i/4e3/2 i28 23/4ei/2 



(3.60) 



Terms that originate from the upper limit of integration in the third interval are expressed in terms of 
half-integer powers gE. These terms (an example is [gEfl'^ = {gEf \/gE involve the square root of 
negative argument because gE < and are thus imaginary (hence, they do not contribute to the real part 
of the result). We have defined the square root to have its branch cut along the negative real axis, and so 
the contributions from the upper and lower branch of the contour C cancel for these terms. Finally, after 
the cancellation of e. To is found as 
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For the second term. 



-g j> dgS2(g) , 



(3.61) 



(3.62) 
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we find 



To 



The third term is 
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The total result of the contour integral of the WKB expansion \sT = Tq + T2 + T4, where 
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The perturbative counterterms read 

P = i ln(27r) - In [F (i - B^iE, g))] + B,{E, g) In (|) 
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where the result [for the definition of B4{E,g), see Eq. (|3.50p ] 

B^{E,g) = B{A, E,g)^E^g(^ + ^E' 
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has been used, as well as the asymptotic expansion 



E' 



1155 
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Inra + z) = z {Hz) - 1} + i ln(27r) - — 



7 



24 z 2880 z3 



(3.63) 
(3.64) 

(3.65) 



(3.66) 



(3.67) 



(3.68) 



(3.69) 



Note that the logarithmic term E hi{—gE/4) in the perturbative counterterms [Eq. p.67p ] is obtained after 
combining the term - lnr(l/2 - S4(£', g)) « E \a{-E) with the term Bi{E,g)\n{g/4) « E ln(g/4). The 
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instanton function A4{E,g) for the quartic potential is finally found to read 

We finally remark tl?at in all intermediate expressions p.52p — p.70p . we l?ave neglected higl?er-order terms 
tl?at contribute in tl?e order g'^ and higlier to tl?e function A4{E,g). The term of order can be found 
below in Eq. (|4~TT|) . 

Alternative methods for the evaluation of the contour integrals about the cut of the WKB expansion are 
described in Appendix F of Ref. [2]. One of these is based on a Mellin transformation with respect to E, after 
which the q integral can be performed with ease, followed by a Mellin backtransformation, which evaluates 
the original integrals. The second alternative method is an adaptation of dimensional regularization and is 
described in Appendix F.7 of Ref. [2]. The method described above, which directly evaluates the integrals 
along contours infinitesimally displaced from the real axis, might be most transparent one among the various 
approaches, and therefore we have chosen it in the current context for illustrative purposes. 



3.5 Generalized Quantization Conditions 



In Sec. 13.11 we have investigated the perturbative expansion for the logarithmic derivative of the wave 
function. The perturbative quantization condition Bm{E, g) = n + ^ was derived for an oscillator of degree 
TO, where m can be either even or odd [see Eqs. p.l9p and p.50p ]. For the "stable" configurations which 
concern even potentials with positive coupling and odd potentials with imaginary coupling Reg = 0, there 
are no nontrivial saddle points of the Euclidean action, and thus no instanton configurations to consider. 
The quantization conditions are thus equal to the perturbative condition, and read 

BM{E,g)^n+^, N even , g>0, (3.71a) 
BMiE,g)^n+^, A/ odd , .g < . (3.71b) 

The question then is how to modify these conditions for those values of the coupling constant where instanton 
configurations become relevant. To this end, let us first observe that we can write the conditions as 

^ =0, (3.72) 



ra-Brr,{E,g)) 



where m can be either even or odd. So, on the one hand, according to Eq. p.47p . we have Bm{E,g) = 
E + 0{g). On the other hand, the quantity l/r{^ — E) can naturally be identified as the spectral determinant 
of the harmonic oscillator with Hamiltonian 



1 92 1 



In this form, we naturally identify 



^--:^%+2«'- (3.73) 
det(-H - E) (3.74) 



2 dq^ 2 

1 



as the spectral determinant. The effect of the instanton is to add an infinitesimal imaginary part to the 
energy of the bound state, 

E ^ 8 = n+\+\lm.E . (3.75) 
Expanding in the nonperturbatively small (in g) imaginary part of the energy, we obtain 

-«-(-l)"+i/2n!Imf (3.76) 
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Now, starting with an even oscillator (for definiteness), the imaginary part of the energy can be approximated 
as [see Eq. (fZ27l) ] 

' ^ ^ exp{-AN{E,g)) , (3.77) 



we have used p.49p . If we now generalize the left-hand side of Eq. p.76p as l/r{^ — BN{E,g)) and 
approximate the right-hand side of Eq. (|3.76p by the right-hand side of Eq. (|3.77p . we have 



1 If 2C{N) 



exp{-Apf{E,g)) 



A^even, 5<0, (3.78) 

which is our conjecture for even oscillators in the "unstable" regime of negative coupling parameter g <0. 
For an odd oscillator, we approximate the right-hand side of Eq. (|3.76p as 

T c T iM)f ^ 1 f 2C{M) f A{M) 



exp (-Am 5)) ■ (3.79) 



2C{M) 



We now generalize the left-hand side of Eq. p.76p as 1//^(| — BM{E,g)) and approximate the right-hand 
side of Eq. (|3.76p by the right-hand side of Eq. (|3.79p . and obtain 

' ^ eM-AM{E,g)], 



r{\-BM{E,g)) V .gi/(*^-2), 

Modd, 5>0, (3.80) 

which is our conjecture for odd oscillators in the "unstable" regime of positive coupling parameter g > 0. 
In writing Eq. p.80p . we correct a typographical error in Eq. (18b) of Ref. [3J, which concerns the missing 

— gi/(M_2) 1 ■ The quantization conditions ()3.78p and p.80p effectively 

assign specific numerical constants to the expression on the right-hand side of Eq. p.43p . 



3.6 Generalized Perturbative Expansions 

We intend to investigate the structure of generalized perturbative expansions which fulfill the quantization 
conditions p.78p and (|3.80p . To this end, we enter into the quantization conditions (|3.78p and p.80p 
with an ansatz that involves an unperturbed energy of the form n + ^ + rj where 77 is a nonperturbatively 
small imaginary part that is assumed to be composed of a nonperturbative exponential factor, multiplied by 
fractional powers of the coupling constant. We then expand the F functions in Eqs. p.78p and (|3.80p about 
their poles and equate the coefficients of the fractional-power terms. In this manner, we easily derive the 
following expressions for the decay width of an arbitrary state of an even or odd anharmonic oscillator. We 
thereby recover the results given here in Eqs. (|2.27p and (|2.29p . but we are able to qualify the order of the 
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relative correction terms. They read 



Imi?W(g<0) 



1 



2C{N) 



1 + 0(52/(^-2) 

^ 2{nN+l)/{N-2) 



exp 



A{N) 



{-g)V(N-2) 



(-g)(2«+i)/(W-2) 
X "1 + 0(52/(^-2)) 
for even oscillators, and for odd oscillators, we obtain 



cxp 



22/(W-2) N_\ 

(_C,)2/(JV-2) 



Ime^)(5>0) 



2C{M) 



I 2[("-l)*^+3l/(A^-2) 
/^n\ g(2n+l)/(2Af-4) 



exp 



A{M) 



9 



l/(M-2) 



exp 



22/(M-2) 



M M 



M-2 ' M-2 ^ 



,l/(M-2) 



1 + 0(5 



l/(J\/-2) 



(3.81a) 



(3.81b) 



Using the dispersion relation (|2.31bp . we recover the result from Eq. (|2.33p for the leading asymptotics for 
the perturbative coefficients of an arbitrary state n of an even oscillator of degree N, 



e: 



(^) (-1)^+1 (jV-2) 

n,K ^3/2 2/^+1-" 
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X r 
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N -2' N~2 
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2 2 



1 + 



(3.82a) 



and clarify that the correction terms are of relative order l/K. The analogous, general result for an arbitrary 
state n of an odd oscillator of degree M reads [see Eq. p.34p . we again clarify the order of the correction 
terms!. 



AM) 



M-2 



7^3/2 ^1 22/f+l-ri 



B 



M 



M 



M-2' AI -2 



-(Af-2) K-n-1/2 



X r 



(M -2)K- 



1 + 



K 



(3.82b) 



The quantization condition ()3.78p implies the following general structure for the generalized non-analytic 
expansion (sometimes called "resurgent expansion") characterizing the resonance energy of the nth level of 
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the A^th anharmonic oscillator (even N), 



J=0 



2C{N) 



exp 



A{N) 



L=0 



2C{N) 



E 



J=0 



(-3)2/(iV- 



2(nN+l)/(N- 
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i-gf 



f:^n\ (-5)(2n+l)/(Ar-2) 



exp 



22/(Ar-2)5( W 



7V-2 ' JV-2 ^ 



(_5)2/(Af-2) 



L=0 



2W/(Af-2) 
(-9)2/(^^-2 



^=0 



(3.83) 



with constant coefficients ^^f^^ ^'^'^ ^ 



= max(0, J— 1). The above triple expansion is complicated and 
in need of an interpretation. The term with J = recovers the basic, perturbative expansion which has only 



{N,n) 

a,0,K{N-2)/2 



(— where the eI^j^ represent the perturbative 



integer powers in g. Therefore, S, 
coefficients from Eq. p.32p . Note that the perturbative contribution to the energy levels is present for both 
positive and negative coupling g. The term with J = 1 recovers the leading contribution in the expansion 
in powers of cxp{~A{N)/g) to the imaginary part of the resonance energy. The term with J = 2 involves 
a logarithm of the form In ^— ^_g^J'/fJ-2) ^ ■ The explicit imaginary parts of the logarithms cancel against 
the imaginary parts of lower-order perturbation series which are summed in complex directions, in a manner 
consistent with the analytic continuation of the logarithm |45l l46l [471 [2]. 



Let us illustrate the expansion ()3.83p by calculating the first terms in the expansion in powers of J up to 
J = 3, including the logarithmic terms, but only up to leading order in g, i.e. only the terms with K ~ 0, 
for the ground state of the quartic oscillator. The result of this calculation is 



Ei^\g < 0) = {i + Oig)] - i exp (^) [1 + 0(,)] 

{1 + Oig Hg))} 



2 / 2 
— exp — 



In 



ng 



3/2 



exp 



■ 7e 



-7e 



■C(2) 



{1 + Oig Hg))} 



(3.84) 



Here, 7e is Euler's constant, and ((2) = -k^ jVl is the zeta function. The terms with J = 1 and J = 3 con- 
tribute to the imaginary part of the resonance energy, whereas the terms with J = 2 adds a nonperturbative 
correction to the real part of the resonance energy which is present only for g < 0. For small coupling, the 
imaginary part is of course dominated by the leading term with J = 1, and the corrections, expressible in 
powers of g, to the (J — l)-term are the phenomenologically most important ones. These correspond to 
the coefficients labelled ^'"iqx ^'^'^ ^''^ given in Eq. (|4.13p below. 

We also note that the quantization condition (|3.80p implies the following structure for the generalized 
non-analytic expansion characterizing the resonance energy of the rtth level of the Afth odd anharmonic 
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oscillator, 
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(3.85) 



K=0 



Here, the perturbation series is recovered as 2p^p']^''^^^_2) = e^'i^ . arid again imax = max(0, J — 1). Just as 
for the even oscillators, it is instructive to illustrate the structure of this expansion by calculating the first 
few terms, restricted to the zeroth order in g, i.e. only the terms with K = 0, for the ground state of the 
cubic oscillator. The result is 
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ng 



exp [l + Oig)] 



— exp 
ng 
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15^ 



In 



- 7E 



155. 
{l + OigHg))} 
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exp - — 
53 



; In 
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{1 + 0{g Hg))} 



(3.86) 



For the perturbative expansion about the first instanton, see Eq. (|4.5p below 
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Chapter 4 

Higher-Order Results 



4.1 Orientation 

In the previous chapters of this article, we have derived the leading-order results for the decay widths 
and the factorial asymptotics of the perturbative coefficients for even and odd anharmonic oscillators. Our 
formulas are applicable to general resonances of anharmonic oscillators of arbitrary degree. In the current 
chapter, we would like to apply the generalized quantization conditions (|3.78p and (|3.80p in order to derive 
concrete results for the higher-order corrections to the leading-order results for the first six anharmonic 
oscillators: those of degrees 3,4, ...,8 (cubic, quartic, quintic, sextic, septic and octic oscillators). The 
material presented in the following sections may seem nothing more than a collection of formulas; yet, 
the evaluation of the higher-order terms is somewhat nontrivial because one has to enter the quantization 
conditions with an appropriate ansatz for the energy in order to derive the higher-order terms. We thus 
felt that it would be useful to indicate higher-order terms for various potentials in order to illustrative the 
applicability of the new results. 



4.2 Cubic Potential 



We 
to 



5 recall that we use the cubic Hamiltonian in the normalization: /i3(.g) = — ^ ^ + 5 <f + \/9'f • according 
Eq. (j2.1b|) . The result for B^{E,g) is easily derived, 

/ 7 15 ,\ , /1365 1155 .\ , , 

B.iE, 9)=E + g[- + -E^)^,^[-^E^^ E^) (4.1) 

3/119119 285285 o 255255 -\ 4 /156165009 ^ 

-I- q 1 E^ ^ E^\ + E 

V 2048 256 128 / V 16384 

121246125 ^, 66927861 .\ ,/ 10775385621 67931778915 o 

H E'^ H E^ \ + q^ [ \ E^ 

2048 1024 / V 262144 65536 

51869092275 . 9704539845 A 

H E^ ^ E^ + O(q^) . 

16384 4096 J ' 

The nonalternating perturbation series for the ground state is obtained as a solution of B3{e'^\g),g) = i: 
(3), , 1 11 465 2 39709 3 19250805 4 5, 
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The higher-order terms of the A function read 

, , 2 /77 141 .\ o /15911 11947 o, , , 

A3{E,g) = — +5 + +5' Ktt^^+^;:^^ (4-3) 



155 " V32 8 y V 128 32 
, /49415863 6724683 . 5481929 .\ , /2072342055 ^ 



V 122880 1024 512 /V 32768 

44826677 . 733569789 A . / 404096853629 
■ ^ ^ KTHZ — ^ ] + 9 



128 2048 /V 1310720 

1100811938289 . 307346388279 . 134713909947 .\ 

H E'^ H E^ H E^^ + O(q^) . 

163840 16384 10240 / ^ 

With C(3) = 4, the conjectured quantization condition for g > 0, see Eq. p.80p . reads for the cubic 
potential, 

1 ^ / R \ ^^(^'ff) 

— — = ^ eM-A,iE,g)). (4.4) 

r(i - B3(E,g)) V87r V .9/ 

The nonperturbative terms can now be determined on the basis of the quantization condition. The solution 
for the ground state, which generalizes formulas found by Alvarez [48J, reads, 

r (3), „^ ^^'Pl^if?) f, 169 44507 2 86071851 o 

Ime^ > 0) = ^ <1 q 

° ' ^/Jfg \ 16 512 ^ 40960 ^ 

_ 189244716209 4 128830328039451 5 1027625748709963623 g 
2621440 41943040 ^ 6710886400 ^ 

933142404651555165943 7 7583898146256325425743381 g 
107374182400 ^ 13743895347200 ^ 

42597573240436340390869984311 „ 

— 5 



1099511627776000 

1509722177 

9^ 



104925715318980163875631509722177 



35184372088832000 

140719287091496360714368165944515013 n i2n1 ^ 

(? +0(0 )>. (4.5) 

562949953421312000 y ^ yy J j v J 

In terms of diagrammatic perturbation theory about the instanton configuration, this result would correspond 
to a 22-loop evaluation of the fluctuations about the instanton (see also Sec. 15.21 belowV This translates 
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into the following corrections to the leading factorial growth of the coefficients, by virtue of Eq. (|2.31ap . 

,(3) ^^J_ f]^Y^^ r(K+^) f 169 64787 48521417 
^■^ TT^/^ \2 J ^ -^2)]^ K 28800^2 5760000 if^ 

418278626369 130056961585097 1459225926418669541 



8294400000^4 331776000000 398131200000000/^6 
5687653597839945011623 66875677760063253018761141 



143327232000000000 K'' 137594142720000000000 K» 
60795686738650254719443177279 



9172942848000000000000 
1906445431654625736765273984543497 



19813556551680000000000000 ifi" 

22298642113252679371826228458403 q ( ^ \\ C4 g-, 

15850845241344000000000 X" ^ [l^ J j ' 

This result confirms, and extends by five orders in the 1/K expansion, the result given in Eq. (11) of Ref. [14J. 
For the first excited state, we have for the decay width up to relative order g^, 

^ (3), , 8c-2/(i5g) r 33349 395368511 , 

Ime] '(g) = ^ <^ 1 g H g'^ 

0F53/2 \ iQ y ^ 40960 ^ 



1788829864593 4 2121533029723423 5 27231734458812207783 g 
2621440 ^ 41943040 ^ 6710886400 ^ ~ 

37583589061337851179291 , 442771791224240926548268373 o 



■9 



107374182400 " 13743895347200 

3458498643307039634422016023691 „ 

9 



1099511627776000 

02225295009 

9^ 



11415793836918299403741902225295009 



35184372088832000 
19846244647223333437533744190161969649 



562949953421312000 



9''+Oig'')}, (4.7) 



where we observe the numerically significant contribution from the next-to-leading order correction term 
— ^ g, even at small g. The first few correction terms have recently been given in [3]. The asymptotics for 
the perturbative coefficients (large K) are 



(3) _ 15^+3/2 3 r 853 135709 426076511 

^i,K - - 2^-3/2^3/2 i +2) I " 120X ^ 28800^2 " 17280000/4:3 

687908048171 3037546625512063 16549184414304538261 



2764800000^4 995328000000 /sT^ 398131200000000/^6 
29911997493767498508937 1417974891656576662560155573 



47775744000000000 K'^ 137594142720000000000 
15158487931720896294670892498891 



82556485632000000000000 

23212260590744971348988221615146323 
6604518850560000000000000 K^^ 



170817014684450099129620368822347752849 1 1 i 

2377626786201600000000000000 /sT" ' ' ' 
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4.3 Quartic Potential 

For the quartic Hamiltonian in the normalization H^i^g) = + +9<f< according to Eq. (j2.1ap . 

the result for B4{E,g) reads 

B,{E, g)=E-g^l + lE^yg' (^^E + ^E^^ (4.9) 

_ /1995 2625 1155 4 

^ V 256 32 16 

. /400785 ^^165]^^3 + 45045 ^,|^^^^,^_ 



V 1024 128 64 

The alternating perturbation series for the ground state is a solution B4_{E)^' (g),g) 



T.(i)r ^ 1 3 21 2 333 3 30885 4 916731 5 6^ ■,r.^ 

The new result concerns higher-order terms of the A function, 

3 / 372101 125333 . 47431 . 
- q-^ 1 E^ H E^ 

^ \ 9216 384 192 

4 /3839943 ^ 82315 . 317629 A 5, 
+ 9 E H E-^ H + O(g^) . 

^ V 2048 16 128 J ^ ' 

The conjectured quantization condition for an even potential at negative coupling g < 0, see Eq. p.78p . 
becomes 



1 1 /4^ Bi{E,g) 



r{\-Bi{E,g)) ^ \g 



eM~ME.9))- (4.12) 



The higher-order corrections to the decay width of the ground-state resonance energy for g < Q can now 
easily be determined, and we again present a 22-loop result (in language of Feynman diagrams, see Sec. l5.2p . 



/ 2 ^^95^ 13259 ,^ 8956043 , 



In.^;-(,<0) = -exp^-j^--^l + -, ^^^2. , ^^g^^ 

11481557783 4 4580883830443 5 12914334973382407 g 
7962624 ^ 191102976 ^ 27518828544 ^ 

6938216714164463905 7 33483882026182043052421 § 
660451885056 ^ 126806761930752 ^ 

201610633027633355940784741 



27390260577042432 



■9' 



296425081490007428403791901173 
1314732507698036736 ^ 

237708867821086878016755601858687 
31553580184752881664 



• (4.13) 
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The corrections to the leading factorial asymptotics of the coefficients thus read, in view of Eq. (j2.31bp . up 
to tenth order, 

4^H-i/^'-N/43-r(i. + i)|i- ''''' '''''''' 



7r3 ^ 27 I 72 X 10368 if2 2239488 if^ 
25875616775 14187855712747 56333049074528335 



644972544 46438023168/^5 20061226008576/^6 
43291321970558012401 301836286153388103809317 



1444408272617472 K'^ 831979165027663872 K» 
2641578170583411485710072165 5664009914699943317155936358237 



539122498937926189056 77633639847061371224064 K^^ 

_ 6633546006727980489585005984230639 (n f ^ \\ 

5589622068988418728132608 /sT" ^ [lO^ ) j ' ^ ' ' 

For the first excited state, the decay width incurred for negative coupling reads, up to relative order g^, 



(4), , / 1 \ / 32 / 371 3371 , 33467903 , 

3^jV"^l' + ^-'~TT52-''+^294^^ 



73699079735 4 44874270156367 5 181465701024056263 g 
7962624 ^ ^ 191102976 ^ 27518828544 ^ 

133606590325852428349 7 850916613482026035123397 o 



660451885056 126806761930752 

9' 



6544779430184823240660057985 g 



9'' 



27390260577042432 

11946916861319310964857723252821 
1314732507698036736 

11601477714905759444960434746130339 
31553580184752881664 



g'' + 0{g^') . (4.15) 



We can now use Eq. (|2.31bp . to infer the corrections to the leading factorial asymptotics of the perturbative 
coefficients, 



371 23341 



72 K 10368 if 2 

36352799 87794214599 67849536512431 356633980356930271 



2239488 is:3 644972544 if^ 46438023168/^5 20061226008576/^6 
350378589485035426237 3041825321871313559291173 



1444408272617472 K'^ 831979165027663872 K» 
32445494418969664996929544513 83262268260174603393926795718797 



539122498937926189056 77633639847061371224064 

114875486821914130002661264320911491 / 1 M 

T^ ■ (416) 



5589622068988418728132608/4:11 \K^^ , , 

For the ground and the first excited state of the cubic and quartic potentials, we have put special emphasis on 
the evaluation of the tenth-order corrections in order to illustrate the potential applications of the methods 
discussed in the current work. 
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4.4 Quintic Potential 

The Hamiltonian is used in the normalization: ft-s^g) = — ^ ^ + ^ 9^ + \/gf': and the result for B^{E,g) 
reads 

o / 118165905 96201105 . 15570555 . 692835 , 

+ E H + H E^ 

\ 8192 2048 512 128 

o / 36358712597025 142306775756145 . 30926063193025 . 

+ 0^ ^ H 

^ V 4194304 1048576 131072 

4140194663605 . 456782651325 . 9704539845 4^ 
■ E H E H A I + 0\g ) 



32768 16384 4096 

The perturbative coefficients describing the ground-state energy diverge rapidly, even in lower orders in g, 
(5), , 1 449 1723225 , 928230645 . 21855598127812155 , 

^" (.9) - 2 - ^ - ^2^ ' 32768 ^ + ^ ' ^' ''^ 

It is of interest to observe that the leading term of the A function can be written in terms of a single 
Gamma function, although from Eq. (|2.16p . one might have assumed that the result necessarily involves 
the expression 5(5/3,5/3) = P^(5/3)/P(10/3). Throughout the remainder of this article, we attempt to 
write the final expressions with the least number of Gamma functions possible, with the Gamma functions 
entering only the numerators. 

Our new result concerns higher-order terms of the A function. Perhaps quite surprisingly. Gamma functions 
enter at orders g^^'^ and g^^^, but not at order g, 

^^./3J!(|^ (^385 ^ 935 ^3 
21/3 73 TT V 32 72 

/21171 132245 , 10865 A 

In view of C(5) = 2^/-^, we have for the quintic potential, 

1 1 / 2^/3 

—j-^ ——— = -=[-—] exp(-A5(£;,.g)). (4.21) 

r [^-B^(E,g)) V8^ V r'^ ) 

The higher-order corrections for the decay width of the ground-state resonance energy read as follows, 

imeo ig>Uj- ^2/3 ^^i/e 7^(25)1/3^ + 2^/^8tt ^ + 

3^6(1) 275^3(1^, 



+ Ui/3 256^2 2i/3 36V3 7r'^ 



ilH^ - ^ + , + 0(//^)') (4.22) 

1152 216 V3 4096 7r3 j ^ ^ > j ^ > 



35 



4.5. SEXTIC POTENTIAL 



The higher-order corrections for the factorial asymptotics of the perturbative coefficients thus read, for the 
ground state, 



28K+2^6K+5/2 V ' 2; ^ UViK 

1 / TT 7r2 275 r^a)\ 1 



^2 



V 84^/3 1176 1176\/3 7r3/ V 504^3 1764 49392^3 

_ ( 1125377 275 \ \ ^ ^ ( ±]] 23) 

U90272 V3 7r3 +49392^2; ^ ) + [k^ ) j ■ ^^'^-^^ 



4.5 Sextic Potential 

With the Hamiltonian in the normalization: HQ{g) — — i ^ + ^Q^ + 9Q^, the result for BQ{E,g) reads 
B6(i?,.g) = £;-.g(^|i? + ^i?3 

2 /21777 5145 . 693 . 

+ E + ^ E- 

^ V 256 32 16 

3 /12746305 ^ 8703695 o 1096095 . 36465 ,^,4, ,^ 

which gives rise to the following perturbation series for the ground state, 

(6), , 1 15 3495 2 1239675 , 3342323355 , 
Ef (9) = -,+^9- ^9^ + 9' / + Oig^) . (4.25) 

Our new result concerns higher-order terms of the A function, 

/2504899 45769 17527 \ 
^ V 7680 96 160 y ^ ^ ^ 

The special role of the sextic potential is illustrated by the fact that by analogy with the quintic potential, 
one would have expected nonvanishing terms of orders g^^^ and g-^/^ in the above result. The surprising lack 
of these coefficients explains the result of a lack of a l/if-correction to the leading asymptotic growth of 
the perturbative coefficients, as is evident from Eq. ()4.29p below. 

With C(6) = \/2, the conjectured quantization condition for the sextic potential reads. 



23/2 

r{^~B,{E,g)) ~ ^ 



6(B,ff) 



exp(-^6(i?,5))- (4.27) 



For the ground-state resonance energy at negative coupling, the following higher-order corrections are easily 
derived. 



Im4''(5<0) = -exp( - 



25/2(_g)l/2; 



1 



2A^/V, 165 174241 2 328657111 ^ 

^ + +^457^^ ^'^(^ )) ■ ^^-2^) 



36 



4.6. SEPTIC POTENTIAL 



Because of the lack of a term of relative order ^/g in this expression, a term of relative order K^^ is thus 
missing from the corrections to the leading factorial growth of the perturbative coefficients. The first few 
correction terms read for the ground state, 



o.A' V ; ^2K+2 V ^2; I 2048 



V 32768 8388608 / \K'^ 



4.6 Septic Potential 

For the septic potential, we encounter the most complex analytic expressions so far in our analysis. The 
Hamiltonian is used in the normalization h7{g) = + ^ + ^/gq"^ ■ The first terms of the result for 

B-j{E,g) read as follows, 

B7iE^g)^E + g('^^'-^E^^'^E^+'^E^ 
^ ^ V 2048 512 128 32 

2 / 182627818702875 156916927352185 . 13513312267455 . 

+ E + E^ H E^ 

^ V 2097152 524288 65536 

824707412529 , 43689020375 9 456326325 .A . 3^ .4 
16384 8192 2048 ^ ^ ■ 

The perturbation series for the ground state is a solution of BT{e'^\g),g) ~ i, 

(7),, 1 44379 715842493569 , 

^" (.9) -2-^^ + ^(^^) • (4.31) 

The new result for the A function contains the golden ratio (\/5 + l)/2, 

A(E a)- I ,1/5 ^2(1)^(1) (^5 9^, 

7^ '^^ 21/10(75+1)1/2 9 7^^1/5 ^ 29/10(75+ 1)1/2 TT 1^8 10 

,3 5Vny5 + i)V2r2(i)r(f) , ^.3\ 

^ 23/10 TT ^1600 2000 y 

+ ,3/5 ^^^^(f)^^(t) /59143 15351 ^ 13209 ^ A _ 

2V10 (75 _ 1)1/2 1^ 9600 400 1000 J ^-J ' ^ ^ 

The quantization condition for g > 0, in view of C(7) = 2^/5, reads 

1 1 



27/5\ BAE,g) 

expi-A,iE,g)). (4.33) 



- 



r{^^Bj{E,g)) 78^ V .g'/V 

The analytic result for the decay width of the ground state reads, in higher order, 

Ime^^Hq > 0) ' exp ( 5V^^(1)/-(|) ^ 

lmeo(g>0)- ^4/5 7^^1/10 I 21/10 9 tt (75 + 1)1/^ gV^ 



\ 1/5 i7r2(|)r(|) / 273(v/5 + i)V2r2(i)r(f) 

- - T 4^ ^ 23/10 53/4 4007r 

+ O (g^/^) \ . (4.34) 



29/1053/4(75 + 1)1/2 \^ 23/10 53/4 4007r 

289r4(|)r2(|) 



24/5 51/2 (75 + 1)1/2 3207r2 
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The higher-order asymptotics to the perturbative coefficients are found as 

5A'+l/2 



22i^+l7r3/2 51/4/^(1) 7^2(2) j 

2^/^17n 1 / ITtt 



51/4(^5 + 1)3/2 225 X is:2 \^ 53/4 21/2 (^+ 1)1/2 900 
17 TT 289 TT^ 

51/4 21/2(75 + 1)1/2 4500 " (75 + 1)405000 



+ 



289^2 9ir4(i)r3(| 



)-(^)} 



(4.35) 



51/2 (75 + 1) 135000 51/2 (75 + 1) 7r2 135000 
are found to be rather involved in their analytic structure. 

4.7 Octic Potential 

As we use the Hamiltonian in the normalization H^{g) = — i ^ + 5 9^+3 9*. we find the result for Bs{E, g) 

^ \128 16 8 

2 / 5604849 ^ 3209745 . 291291 6435 ,\ 3, 

+ g^ E-\ H £'-^H E^]+0(g'^). (4.36) 

V 2048 512 128 32 / ^ ^ 

Solving Bs{EQ^\g), g) — ^, we find for the ground state: 

^(8), , 1 105 67515 2 401548875 , 25424096867715 4 ^ , 

4 ^ (.9) = 2 + 5 - ^ .9^ + 5^ ^^5^^ 9' + Oig-) . (4.37) 

The result for the subleading terms of the A function is much less involved as compared to the septic 
oscillator and reads, 



^8(^'ff)-22/3io,(„5)i/3+(--9) l^n^ (,16 + 4^ 

-T^d) (n , 91£:2\ ^28007 , 22669 ^2 , 2587, 



The coefficients of the terms of orders g^^^ and g^/"^ involve Gamma functions, whereas the term of order g 
is a polynomial in E with rational coefficients, in analogy to the corresponding result for the quintic oscillator 
recorded in Eq. (|4.20p . In view of C(8) = 2i/'', we have the following quantization condition for 5 < 0, 



1 1 



/24/3\ 



B8(E,g) 



eM-ME,9)). (4.39) 



r{^~Bs{E,g)) 72^ 
By virtue of the dispersion relation (|2.31bp . the following result for the decay width in higher order, 

73r3(i) \ 1 / 2^^/^ 



(8)^„/n^-_... ^3^ 



ImE^^'ig < 0) = -exp 



22/3l0 7r(-g)i/3 J \ g, 

2\ / /lorSi'lA 101 r6/'2^ 



1 - (-,)i/3 ^^^'(i) + (-,)2/3 f^n^ + ^^(i) ^ 

^ 21/3 8 73 TT ^ ^2/3 108 73 ^ 22/3 384 712; 

, , /12429 5397r 133ir9(|) \ 

' ' ' ■ ■ ) +0{g^/^)\ , (4.40) 
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translates into the following corrections to the leading factorial growth of the coefficients: 

3K+1/2 



o.if V ; ^3/2 \ril)y ^ ' ^' r 180 VSK 

V 1080 V5 194400 ^ 129600 VSttS/ ^ {k^ > ^ ^ ' ' 



39 



4.7. OCTIC POTENTIAL 



Table 4.1: General properties of even and odd anharmonic oscillators of arbitrary degree (even integer N 
and odd integer AI). We summarize the most important interconnections between different regions in the 
complex coupling parameter plane. 



Odd, Unstable 


Even, Unstable 


Odd, Vr (Stable) 


Even, Stable 


(positive coupling) 


(negative coupling) 


(negative coupling) 


(positive coupling) 
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Hn(-9) 


''"'^-'^ = -2 0^ 


1 




1 2 1 1 iV 


1 2 I ., ^1 M 

+ -<? ±i\V9\q 


+ '-^ + 191," 


Under complex scaling: 


Under complex scaling: 


No resonances 


No resonances 


Resonances appear 


Resonances appear 






Perturbation theory is 


Perturbation theory is 


Perturbation theory is 


Perturbation theory is 


not Borel summable 


not Borel summable 


Borel-Leroy summable 


Borel-Leroy summable 


in the ordinary sense, 


in the ordinary sense, 


to the real eigenvalue 


to the real eigenvalue 


but Borel— Leroy summable 


but Borel-Leroy summable 






in the distributional sense 


in the distributional sense 
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Chapter 5 

Views from Different Perspectives 



5.1 Even and Odd Oscillators: An Overview 

In order to view the investigations reported here from a broader perspective, let us remember here that 
originally, the investigations on anharmonic oscillators started with the consideration of the cubic and quartic 
cases as paradigmatic examples for the lowest-order odd and even anharmonic oscillators. One can the 
generalize these in many ways: to an internal symmetry 0{N) [38], to more space-time dimensions (toward 
field theory, see [49,^0]), and to anharmonic terms of arbitrary degree [20,£]. The latter generalization is 
the one studied by us here. In the broader context of the various generalizations mentioned, the dispersion 
relations (j2.31ap and (|2.31bp might be of particular interest [14]. It has often been asked whether cubic 
theories and those with higher-order odd perturbations represent physically self-consistent theoretical models. 
Finally, their Hamiltonians are not bounded from below for field configurations which accumulate near the 
region where the potential assumes large negative values. Based on the considerations presented here, we 
can say that if one continues the coupling into the complex plane, the "odd" theories are no less natural 
than the "even" ones: namely, the "stable" regions for the coupling parameter g (the domain of positive g 
and the domain of purely imaginary g, respectively) are connected to the "unstable" regions (negative g for 
the "even" theories and positive g for the "odd" theories) via dispersion relations. 

One can ask the question why we construct the dispersion relations (|2.31ap and (|2.31bp with resonances, 
i.e. by assigning to the "unstable" parameter regions resonance eigenvalues with a nonvanishing imaginary 
part. The corresponding eigenfunctions fulfill complex-rotated boundary conditions |5]. Finally, one can also 
associate L^(]R)-boundary conditions to the quartic potential, which leads to a spectrum that is unbounded 
from below but discrete (see [51] for an instructive discussion). The reason is that every real energy eigenvalue 
has a natural analytic continuation to the complex coupling plane only when the boundary conditions are 
rotated, i.e. in that case it is possible to find a spectrum of resonances whose real parts are in fact bounded 
from below, even for the "unstable" parameter regions. As stressed in [14], the dispersion relations associate 
the Herglotz properties of the resonance eigenvalues with some global analytic properties in the cut complex 
coupling plane. 

The most important properties of the anharmonic oscillators studied here are summarized in Table 14.11 
A few literature references regarding the entries are in order. Recent numerical investigations regarding 
the resonances in the cubic and quartic potentials, for the the unstable parameter regions, can be found 
in [52] EH [54], and important steps toward a strong-coupling analysis of the resonance eigenvalues of the 
cubic oscillator have been published in [55^ |56J. For the cubic, T^T-symmetric Hamiltonian, the boundedness 
of the spectrum (from below) and the discreteness have been proven in [571 [58]. According to [17], there 
is numerical evidence that the perturbation series for the cubic, T'T-symmetric Hamiltonian is Stieltjes, and 
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Figure 5.1: Comparison of the analytic results (|5.ip and (j5.2p with numerical 
data for the imaginary part lme^^\g) of the resonance energy of the first 
excited state of the cubic potential. 



that the results of Pade prediction techniques are indistinguishable from those for the quartic oscillator, 
which is known to be Borel-summable. Note that the Stieltjes property alone would already imply Borel 
summability. We would like to conjecture here a generalization of this property to all the odd, PT-symmetric 
anharmonic oscillators: namely, that their perturbation series are Borel summable to the energy eigenvalue. 
This conjecture is inspired by the fact that recently [59J, the distributional Borel summability [60] has been 
proven for odd, unstable anharmonic oscillator of an arbitrary degree. We also recall that the dispersion 
relation joining the odd, unstable oscillators with the odd, stable oscillators has been given in [T4], whereas 
the corresponding relation for the even oscillators has been derived much earlier [39]|20j[6]. We also recall 
the proof of the Borel summability of the perturbation series for the quartic oscillator [3- The necessary 
generalization of Borel summability to Borel-Leroy summability for oscillators of higher degree has been 
discussed recently [61j. 



5.2 Numerical Calculations for Weak Coupling 

Let us note first of all that all analytic results derived here have been subjected to extensive numerical 
verification. Here, we would only like to discuss one particular case, the higher-order formula for the subleading 
corrections to the decay width of the first excited level of the cubic potential, which we recall here from 
Eq. (EID, 

T (3)/ N 8e-2/(i59) r 853 , 33349 2 395368511 3 , 4,! ,^ ^, 

Our purpose is to compare this result to the numerical values generated exclusively by the leading-order term 

8g-2/(15g) 



Im 4jcading(5) = - °^/^ 3/2 ■ (5.2) 



The relative correction of order g as given by the term — 853.g/16 « — 53g in Eq. (|5.ip is extremely important 
in reaching satisfactory agreement of the analytic result and numerical data for the imaginary part of the 
resonance energy lme'f\g), even at relatively small coupling parameter g < 0.01. This is evident from 
Fig. 15.11 where we compare the leading term of the asymptotics given in Eq. ()5.2p as well as the first three 
correction terms given in Eq. (|5.ip to numerical data. 
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5.3 Numerical Calculations for Large Coupling 

As we look at the dispersion relations (|2.31bp and (|2.31ap . we might ask ourselves how far the cut extends: 
from g = to 5 = oo, or possibly up to a finite value of the coupling. Our conjecture here is that the 
resonances of the unstable oscillators persist for any value of the coupling, however large. Intuition might 
suggest otherwise. The unstable case is naturally associated with tunneling through the potential barrier. As 
the coupling increases, the barrier becomes narrower until, in the limit of very large modulus of the coupling, 
it eventually disappears. One might assume that only a limited number of resonances persist for a large, but 
finite value of the coupling: namely, only those for which the real part of the energy lies significantly "under 
the top" of the barrier. Numerical evidence [54] suggests otherwise. Namely, we investigate the strong- 
coupling expansion for the resonances of odd oscillators. After a scaling transformation q — >• 5-1/(2(^^+2)) 
which leaves the spectrum invariant, the odd Hamiltonian (j2.1bp reads 

hsMg) = 9'"^"^'^ {hLM + ^ 5-^/'*'+^^ q') , (5-3) 

where hi^M does not depend on the coupling g. This Hamiltonian now gives the strong-coupling expansion 
in a natural way: One can regard the term i 5^2/(*^+2) as a perturbation of the term q^^ . Via a complex 
scaling transformation g — > q&^^\ we transform h£,M — hi,M{0) where 

A diagonalization of this complex-scaled Hamiltonian in the space of harmonic oscillator eigenfunctions 
gives us approximations to the resonance energies. The resonance eigenvalues have a well-defined complex 
argument, and they can be numbered by an index n — 0,1,2, They read 



(M) _ (M) -i7r/(M+2) 



where by virtue of (|2.31ap . we associate the resonance as opposed to the antiresonance with the region 
infinitesimally above the real axis, 



^^(5 + iO) - [eZ^ + O (5-2/(^^+2))] , (5. 



for 5 — > 00. These asymptotics suggest that all resonances n = 0, 1, 2, . . . persist for any modulus of the 
coupling. Specifically, according to Table 2 of [54J, we have 

efl = 0.762 851 775 e^''^/^ , tfl = 2.711 079 923 e"'''/^ , ef^ = 4.989 240 088 e"'''/^ . (5.7) 

The above results define the large-coupling coupling asymptotics for the for the first three resonances 
^n\g) — >■ 9^^^ (-^tn' ^o'' 9 ^ 00 + iQ (with n = 0,1,2). We have also obtained corresponding results for 
higher excited states and conjecture here that all resonances of the cubic oscillator persist for any value 
of the coupling, however large. The corresponding large-coupling asymptotics for the quantic potential are 
given by the levels 

e(^|) = 0.709 935 516 e-'''/^ e^^i' = 2.659 756 364 e-'^/^ = 5.458 235 423 e"'''/'^ . (5.8) 



For the even oscillators, the corresponding formulas read as follows. We recall that the resonances appear 
for negative g, and that the resonance as opposed to the antiresonance is attached to the negative g axis 
for infinitesimal positive imaginary part of the coupling. Under a scaling transformation q (—.9)^^^'^^^^ q. 
the Hamiltonian (j2.1ap transforms to the following scaled Hamiltonian with the same spectrum. 
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As a WKB-approximation to the wave function shows, the appropriate scaling transformation for even N is 
also of the form g — > q&^'^ , but with He.Nid) and On given by 



1 ^ 

2 dq^ 



27r 

iV + 2 



(5.10) 



The eigenvalues again have a well-defined complex phase, 



E. 



(N) 



E 



(N) 



-2i7r/(Af+2) 



(5.11) 



The large-coupling asymptotics for the resonances then reads 



i?W(ff + iO)^(-5)^/(^+^) 



E. 



(N) 



O g 



-i/iN+2)\ 



(5.12) 



for g — >• — oo. A numerical determination of the complex resonance eigenvalues E^^^J of the non-Hermitian 
Hamiltonian He^N{0) is thus sufficient in order determine the large-coupling asymptotics of the resonance 
eigenvalues of iJs.A'(g), which coincide with the resonance eigenvalues E^\g + iO) of H^lg) attained for 
g <0. For the first three resonances of the quartic potential, we find by numerical calculations. 



= 0.667 986 260 e-'''/^ 



E^^l = 2.393 644 016 e~'^/^ , 



Ej,^^ = 4.696 795 386 e"'''/^ 



(5.13) 



This example, which demonstrates the persistence of the quartic resonances in the strong-coupling limit, 
concludes our analysis of the complex resonances. Note that the strong-coupling expansions, whose lead- 
ing terms are indicated in Eqs. ()5.6p and ()5.12p . are in some sense complementary to the weak-coupling 
expansions (|3.83p and (|3.85p . 
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Chapter 6 

Conclusions 



Anharmonic oscillators form a very basic model for many physical processes. The anharmonic terms induce 
coupling among the energy levels of the harmonic oscillator, which can be interpreted as inducing some 
sort of "higher-harmonic generation" due to the intertwining of an infinite number of harmonic modes 
in the energy levels of the anharmonic oscillator. At the same time, the anharmonic oscillators represent 
classically integrable systems. There is no chaos, and therefore, the corresponding quantum systems are 
amenable to a unified treatment which is based on semiclassical expansions and which is manifest in the 
generalized quantization conditions (|3.78p and (|3.80p . These quantization conditions allow us to determine 
the subleading corrections to the decay widths of arbitrary levels of anharmonic oscillators of arbitrary degree. 

The formulas for the subleading corrections to the decay width of the resonances, as presented in 
Eqs. (|4T1) . (|4a3l) . HjAAb^i . (g^ll). (14^281) . (|434)) and (14^401) have not yet appeared in the literature to 
the best of our knowledge. We note, in particular, that the results for the excited states in Eqs. (|4.7p and 
()4.15p are derived on the basis of the generalized quantization conditions; these could not have been ob- 
tained by standard field-theoretical investigations which involve a higher-order perturbative expansion about 
the instanton configurations, as it is the case for the ground-state results (see also Sec. 15. 2p . All analytic 
results presented here have been subjected to intensive numerical checks, including those for the sextic 
anharmonic oscillator where the higher-order corrections are observed to have a rather surprising structure. 

Let us now attempt to interpret the results from a broader perspective. In the case of a theory with an 
"unstable" potential where, classically, the motion of the particle is not confined to a finite spatial region, 
one typically encounters quantum mechanical resonances with a nonvanishing decay width, when one endows 
the problem with "outgoing" boundary conditions (and the corresponding antiresonances for the complex 
conjugated boundary conditions). The complex boundary conditions lead to a non-Hermitian Hamilton oper- 
ator whose resonance energies are not real. Intuition suggests that there must be way of smoothly changing 
the coupling constant and the boundary conditions from the "resonance case" to the corresponding complex 
conjugated antiresonance, so that the resonance eigenvalue "transits" through the real axis. At this point, 
precisely, at least for the odd anharmonic oscillators, the eigenenergies become real, and the Hamiltonian 
becomes 7'T-symmetric (see [15]). This is the observation that leads, naturally, to the dispersion relation 
()2.31ap for odd anharmonic oscillators. 

Analogously [39][2Ql[6], for even anharmonic oscillators, if one starts from negative coupling, there is a way 
of changing the complex argument of the coupling parameter so that the boundary conditions for resonances 
and antiresonances deform smoothly into each other, crossing an entirely real eigenvalue for the energy at 
positive, and real, coupling g, and this corresponds to the relation (j2.31bp . 

In some sense, T'T-symmetric odd anharmonic oscillators receive an interpretation as the analogy complete 
analogy of even anharmonic oscillators for positive and real coupling. They provide a bridge between reso- 
nances and antiresonances, as exemplified by the dispersion relation (j2.31ap . over the real axis. They also 
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facilitate the interpretation of tlieories wfiich allow for the presence of complex resonance and antiresonance 
energies, by allowing us to interpret these theories as the unstable realizations of a stable, 7'T-symmetric 
theory, continued into the complex plane. 

By a similarity transformation [151 [21] [23l [24] , it is possible to map the dynamics induced by a T'T-symmetric 
theory onto a Hermitian theory, but the resulting expressions can be very complicated, and even in the case 
of the relatively simple cubic PT-symmetric oscillator, the corresponding Hermitian Hamiltonian obtained 
by the similarity transformation cannot be written down in closed form and contains operators that mix 
momenta and coordinates. In that case, the original, T'T-symmetric theory, namely the T'T-symmetric 
cubic oscillator, provides for a much more compact formulation that allows for a workable, and practically 
feasible formulation in terms of variational principles and corresponding classical trajectories of the dynamical 
variables, which would be much more difficult to achieve when starting the discussion from the much more 
involved expressions that occur in a typical case [24J for the equivalent, Hermitian theory. This is one of the 
reasons why PT-symmetric theories have received some attention in the physics community in recent years. 

Although a high degree of unification can been achieved for odd versus even anharmonic oscillators using 
concepts from "PT-symmetry and generalized quantization conditions, let us finally, that there is some 
room left for future investigations. E.g., one might well ask what the quantization conditions (|2.31ap 
and (|2.31bp should look like for arbitrary complex argument of the coupling g. Our conjectures as pre- 
sented in Eqs. p.78p . (j3.80p and p.7ip are relevant only for either positive or negative, or purely imaginary 
coupling parameter. These parameter ranges represent the physically most important scenarios by far, but 
we would like to point out that the interpolation required to describe an arbitrary complex argument of the 
coupling parameter remains an open problem. 

Finally, in view of the phenomenological significance of anharmonic oscillators, we believe that the numerically 
large correction terms in the higher-order formulas for the decay widths, as obtained, e.g., in Eqs. (|4.7p and 
(|4.15p . might find interesting applications in the various branches of the theory of finite quantum systems. 
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